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We calculate the exact autocorrelation exponent A and persistence exponent 9, and also ampli- 
tudes, in the dilute limit of phase ordering for dimensions d > 2. In the Lifshitz-Slyozov- Wagner 
limit of conserved order parameter dynamics we find 9 = 7^6, a universal constant times the volume 
fraction. For autocorrelations, A = d at intermediate times, with a late time crossover to A > d/2 + 2. 
We also derive A and 9 for globally conserved dynamics and relate these to the q —> 00-state Potts 
model and soap froths, proposing new poisoning exponents. 
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While much has been learned about the coarsening ki- 
netics that follows a temperature quench from a single- 
to a multi-phase state [Q, relatively little has been es- 
tablished for certain recently introduced dynamical expo- 
nents. It is accepted that the characteristic length scale 
of strongly correlated regions grows as a power law in 
time, L ^ t^^^, with universal z [||. Most non-conserved 
order parameter systems yield z — 2 and those with 
scalar conserved order parameter dynamics yield z — 3, 
independent of the system dimensionality d or of con- 
served quantities such as the volume fraction e of the 
minority phase. Consequently, persistence ^-p^ and au- 
tocorrelation ||l^[r^ exponents, 9 and A respectively, are 
being explored in the hope that they contribute to a char- 
acteristic set of universal exponents analogous to those of 
equilibrium criticality. However, fundamental questions 
remain about the universality of these new exponents, 
and even of the existence of power laws in the relevant 
quantities. 

The persistence exponent 9, introduced in the exper- 
imental study of breath figures [Q, is defined by the 
power-law decay, P{ti,t2) ^ ^2^^, of the probability that 
a stochastic variable has not crossed some threshold — 
typically its mean — between the times ti and ^2 (D . We 
consider the persistence of a local, scalar order parameter 
(/)(r, t) (rescaled so that in equilibrium (j) = ±1), given by 
the fraction of the system that has not undergone phase 
change between ti and ^2 0|- 

In d = 1, 9 has been calculated exactly for the non- 
conserved q-state Potts model ||^, and has been shown 
to be universal by renormalization group methods in the 
Ising {q — 2) case In higher dimensions, studies have 
focussed on diffusion models Q , which exhibit non-trivial 
values of 9, and on related Gaussian approximations for 
non-conserved ordering kinetics [@j8|. For d = 2, these 
approximate results compare well with simulations PpO| 
and with experiments on twisted nematic liquid crystals 
[ pl| , an Ising analog. However, there have been no pre- 
vious studies of persistence for conserved coarsening dy- 
namics. 

In phase-ordering systems the autocorrelation func- 



tion, A{ti,t2) = {(I){r,ti)(f>{r,t2)) — (0)^, decays asymp- 
totically as A(t\,t2) ^ \L(t\) I L(t-2)'^ , which defines A 
||l|,^. For non-conserved scalar coarsening, A has been 
measured experimentally in d = 2 jl^ and calculated 
in d = 1 Ip^ . Approximate calculations and numerical 
results have been obtained for various e and d [T5| ] in 
the case of globally conserved dynamics — i.e., a non- 
conserved order parameter (hence z = 2) subject to a 
constraining field that maintains the total volume frac- 
tion of each phase. Locally conserved dynamics has been 
studied numerically in two (l^Jl^ and three ||l^ dimen- 
sions, and a formal asymptotic bound, A > d/2 + 2, has 
been established by Yeung, Rao and Desai (YRD) [p^ , 
but otherwise little theoretical progress has been made. 

In this Letter, we study persistence and autocorrela- 
tions for both locally and globally conserved dynamics. 
We focus on the asymptotic, late stage regime which fol- 
lows a quench to a subcritical temperature T < Tc, in 
the limit of vanishing volume fraction, e — > 0^. This is 
the limit of the classic Lifshitz-Slyozov- Wagner (LSW) 
theory, which firmly established z = 3 for dilute, locally 
conserved coarsening [p]Jl8|. We use LSW theory in a 
similar spirit, to obtain 9 and A. A summary of our re- 
sults follows. 

For locally conserved dynamics we compute 9 for all 
d > 2, and demonstrate that (i) the persistence decays 
as a power law, P{ti,t2) ^ (^1/^2)^, (h) the exponent is 
a function of the volume fraction, going as 

9 = ide (1) 

in the small e limit, and (iii) is universal in that it does 
not depend on the surface tension, quench depth, temper- 
ature, or mobility, with 72 ~ 0.39008, and 73 ~ 0.50945. 
A large-c? expansion gives the asymptotic series 

k 

7d = y3d78^[l -I- a^d-"' + 0(d-'=-i)], (2) 

m— 1 

which is quite accurate in d = 2, 3 when truncated at 
fc = 3 ||l^jp9. We also compute the order e'^/^ correc- 
tions to dlD in d = i. 
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For the autocorrelation function we find A = d as 
e ^ 0^, with exphcit, universal expressions for both the 
amplitude and logarithmic corrections. We also present 
a physical scaling argument that predicts a crossover to 
X > d/2 + 2 for any finite after t2 > 6'^/%, thus 
satisfying the YRD bound [|l6| . 

Next, we consider globally conserved (GC) dynamics, 
again in the small e limit, and find 9 to have the same 
form (0), with a different, universal 7^*^- In particular, 
0.48797 and 7^ ^ ~ 0.62450. In large d 



72^^ 



ld-' + §-sd-' + 0{d-% (3) 



which is highly accurate for all d . The large d asymp- 
tote, 9 ~ OAOeVd is similar to the approximate result 
9 ~ O.lbVd obtained at e | (where the GC dynam- 
ics is equivalent to non-conserved). The autocorrelation 
exponent is A = c? , with no crossover expected at late 
times. We also find universal amplitudes and logarithmic 
corrections. 

Finally, we draw a connection between the above re- 
sults for persistence in GC dynamics with poisoning (de- 
fined below) in soap froths. First, there is some evidence, 
mainly in d = 2, that soap froths have the same asymp- 
totic dynamics as the non-conserved q-state Potts model 
in the q ^ 00 limit Second, the coarsening of the q- 
state Potts model and that of GC dynamics with e = 1/q 
were shown to be equivalent as g — *■ 00, within a Gaus- 
sian approximation scheme p^ ]. However, the details of 
the topological rearrangements may be different between 
Potts models and soap froths ||2l[], and further, while both 
Potts models and soap- froths have vertices, GC systems 
do not. 

Nevertheless, numerical studies found little difference 
in A between the Potts and GC models ||l5|, implying 
that they may lie in the same dynamic universality class. 
Also, the Potts persistence exponent, which measures the 
volume never visited by a wall, is given via the GC cor- 
respondence hy 9 = d/2 asg— s-cx) [ p2| , consistent with 
Potts simulations [pp^ and with experiments on d = 2 
soap froths |2^. To further explore these analogies, we 
define a new set of poisoning exponents, 9^, that give 
the decay of volume that has never been visited by any 
of a set of phases that occupy a total volume fraction 
S. This poisoning should provide a more delicate test of 
the underlying dynamics than autocorrelations or persis- 
tence, and may be directly explored via simulations of 
Potts models and experiments on foams |2^]. By use of 
the GC correspondence we obtain, for S <C 1, 



7d ^■ 



(4) 



We begin with LSW theory, which applies to widely 
separated drops of minority phase as e ^ 0+. This the- 
ory provides the only solution of a phase-ordering sys- 
tem in d > 1 with topological defects — in this case 
with domain walls. Drops of radius R evolve according 
to R^R — ad{R/ Rc^ 1), where ad is a nonuniversal con- 
stant [0, and the dot indicates a time-derivative. Here 



Rc = i^Oidty^^ ^ L is the critical radius, where drops 
shrink for R < R^ and grow for R > R^. The density 
n{R, t) of droplets of size R at time t obeys a continuity 
equation h = —dfi{Rn), which leads to a scaling solution 
UdiR.t) = R,{t)-d-^fdiR/Rc), with 



eFdX^ exp[-d/(3 - 2x)] 
(3 + a;)i+4'i/9(3-2x)2+5'i/9^ 



(5) 



and fd = for x > 3/2 The normahza- 

tion constant Fd is determined by the volume fraction 
of the minority phase, e = Jq dxfd{x)x'^, where 
Vd = 7r''/2/r(l -I- id) is the unit d-sphere volume. This 
gives F2 ~ 37.752 and F^ ~ 186.13, and the large-d ex- 
pansion VdFd = 6^^21/2+8^/9 y27d7^[H-S^^i a™d-™ -|- 
0(d~'^~i)] [^. The total number density of drops can 
be shown to be n{t) = eF<j/(4d3'^e'^/3)[i?c(t)]-'*. 

The droplet growth equation can be written in terms of 
the scaled size x = R/Rc{t) as itx^x — —{x + 'i){x — |)2; 



where i < for all x > 0. 
give the trajectory 



This may be integrated to 



t{l 



(1 



fx)5/3exp[2x/(3-2x)], (6) 



which is parameterized by the time of complete evapora- 
tion tf. From (^, we find an expansion for ti t2 < tj 
which will prove useful below: 



X{ti) = I [1 - S{ti/t2) + 



(7) 



where 5(^1/^2) = 1/ [ln(t2Ai) + f lnln(i2/ti)] . This 
leading correction to x{ti), related to the essential singu- 
larity in fd{x) at X = |, is universal and independent of 

X{t2). 

Using these LSW results, we can calculate P<(ti,t2), 
the persistent or unpoisoned volume fraction of minority 
phase, and P^{ti,t2), that of the majority phase. The 
total persistent volume P = -\- will decay with 
the slower of the two unpoisoned fractions. 

Only droplets that have survived to time t2 contribute 
to P<, the unpoisoned minority phase, and their den- 
sity decays as n ~ R^*^ ^ ^2'^^^ ^ using the notation 
Ri = Rc{ti). The density of droplets that have survived 
and are smaller at t2 than they were at ti is 



Ro 



dxfdix) 



R-'^-^ ^ t-"'-\ (8) 



where we have used the small x behavior fd{x) ~ x'^ for 
t2^ ti, and Rc ^ t^^^. Hence, droplets that have shrunk 
comprise a vanishing fraction of those surviving at ^2, so 
that no surviving droplets have been poisoned at time 
t2 ^ ti. Consequently, P< is asymptotically the num- 
ber density n(t2) times the volume at the initial time ti 
of those largest droplets, Vd[x{ti)Ri]'^. Using x{ti) from 



P<(ti, i2) = iePd(tiA2)''/'[l - d6{t,/t2) + ...], (9) 
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where Bd = FdVa/ with B2 ~ 7.6115 and 
B3 ~ 11.951, and the leading logarithmic corrections 
5{ti/t2) are universal. 

The majority phase, with unpoisoned volume fraction 
P^, can only be poisoned by growing drops, i.e. those 
with R > Re- Since the drop positions are uncorrelated 
in the dilute limit ||l^ — the key feature which makes 
LSW theory soluble — it follows that the unpoisoned 
regions must be uncorrelated as well, leading to 



(10) 



where v{t) is the rate of encroachment by minority phase. 
From LSW theory we have 



where the second equality comes from mass conservation 
of drops larger than Rc{t), i.e. dt J^^^ dx x"^ fd{x) = 0. 



Using Rc oc i^/^, we obtain v 
where 



Vdfdil)/m = ejd/t 



Id = FdYd/iA'+^^^'Se") 



(12) 



is a universal constant. Combining ( p^ and the initial 
condition, P^(ii, ii) = 1 — e, we find 



F>(ti,<2) = (l-e) iti/hr'\ 



(13) 



so indeed decays as a power law. Remarkably, this 
result is valid for all t2 > ti in the scaling regime, not 
just when ^2 ^ ^i- As expected, P-^ decays slower than 
P< , and so P(ti,t2) ^ P^{ti,t2), leading to equation 

In order to derive A from the persistence, it is conve- 
nient to change field variables to i/j — (0 + l)/2 (with 
minority phase ip — \ and majority ■(/; = 0), giving 
(-0) — e and the autocorrelation function A{ti,t2) = 
4[(V'(r, ii)-i/;(r, t2)) — e^]. The two-time average is then 
the probability of finding a given point inside minority 
droplets at both t\ and t2- The contribution from un- 
poisoned regions is exactly P<(ii,t2), whereas poisoned 
regions that find themselves in a droplet again at t2 con- 
tribute [e — P^{ti, 12)]'^. To leading order in e, 



A{h,t2) = eBd{Ri/R2f[l - dd{h/t2) 



(14) 



giving A d in the dilute limit (this was noted before in 
the GC case fl^). This exponent depends solely on the 
existence of a scaling distribution of uncorrelated drops; 
in the LSW limit the details of the drop distribution and 
evolution serve only to determine the universal amplitude 
and leading logarithmic corrections. 

We expect (^^ to hold for R2 much less than the drop 
separation at ti, Ri ^ e~^^'^Ri. For e ^ 0+ this is for- 
ever. Below we discuss e > 0, where the drop separation 
is finite and correlations must be considered. 



dR VdR'^dtn{R, t) ^ VdRin{Rc, t)Rc, (11) e^'^ ^fuj^ [l 



We turn now to globally conserved dynamics, where 
the analog of LSW theory was given by Sire and Majum- 
dar 115). Droplet growth follows R = ad[l/-Rc(i) - 
where Rc{t) represents a time-dependent applied field 
tuned to maintain the volume fraction, e. Combining 
the droplet growth with the continuity equation, as in 
LSW, we find a scaling solution for the droplet density, 
niR,t) = R-''-^f(^^{R/Rc), when R^it) = {^adtY^^ 
with the distribution 

/GC(2.) ^ eF2^x(2 - x)-'^-2 exp[-2rf/(2 - x)] (15) 

for x < 2, and fd ^ otherwise. The normalization 
condition, e — Vd dx x'^ f^'" {x), determines P^*", 
with Pf ^ ~ 16.961 and F^^ ~ 120.29. A large-d ex- 
pansion yields the excellent approximation V^P^*^ = 

d-'+§,d--~mi,d-'+o{d-% 

The droplet density is {t) = ePf^ {2e)~'^ / d [|o|. In- 
tegrating the scaled growth equation, 2txx — — {2 — x)^, 
gives the trajectory tf — t{l — ^x)^ exp[2x/(2 — x)]. 
Hence, drops surviving to time t2 ^ have 



x{ti) = 2[l-5'-'-{ti/t2) + 



(16) 



where the leading correction is 6^^{ti/t2) — 
2/ [ln{t2/ti) + 21nln(i2/ii)]. The density of drops that 
are smaller at t2 than at ti decays as i?^''^^, so these 
are again negligible asymptotically. Consequently, the 
autocorrelation function A{ti,t2) — 4P<(ti,t2) is 



AihA2) - ei?^^(i?i/i?2)''[l - dd''^ih/t2) + 



(17) 



with B^^ = VdF^^/{de'^), giving B§^ ~ 3.6057, 
B^'-' ~ 8.3623. The leading logarithmic corrections in 
5'^^[ti/t2) are universal. 

The calculation of majority poisoning and persistence 
goes through the same as before, using (^ and (p^), 
with the result 9 = l^'^e. The different growth exponent 
i?e - t^'^ gives 7d^^ = yd/f^(l)/2e = ^^^^^7^-2^/2. 
The result is (|). 

Up to this point, calculations of, say, /^(a;), have been 
for the leading 0(e) term, for which the drops may be 
regarded as uncorrelated. Higher order effects such as 
droplet collisions and diffusion-mediated interactions will 
lead to correlations in the drop sizes and positions. How- 
ever, screening of the diffusion field has been shown in 
d = 3 to contribute 0{e^/^) corrections for uncorrelated 
drops [p7|-p9|, which is thus believed to represent the 
leading correction to LSW theory, with correlations com- 
ing in only at O(e^) Note that with GC dynamics 
there is no diffusion field, hence the leading corrections, 
due to collisions, are expected to be of O(e^). 

With our existing machinery, then, we can compute 
the leading corrections to the LSW exponents. Since 
drops are uncorrelated to 0(e'^/^), \ — d remains un- 
changed. However, the distribution function is f3{x) = 
h{x)[l -f ei/2{G'3 + 53(x)}] + 0(e2), with 



53(2:) = bo 



2 In 



64a; - 87 
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+ 3 4(3-2a;)2 



(18) 
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where 60 = IjAWe ~ 1.6297, and we have maintained 
X = R/Rc P0|] . The normahzation condition determines 
G3 ~ —3.4047. This leads, via ( pTj) , to the persistence 
exponent 9 ~ 0.50945e - 0.14969e^ + 0{e'^). 

The leading correction to A for e > is of a different 
nature. At sufficiently late times t2, A becomes strongly 
dependent on correlation effects [Q. This occurs when 
drops grow to be larger than their earlier spacing, R2 ^ 
Ri ^ e^^/'^Ri. Then, the autocorrelation function is no 
longer given by the unpoisoned minority volume fraction, 
since each drop at t2 covers many drops from ti . As a re- 
sult, the decay of autocorrelations becomes dominated by 
the fluctuations in drop density at ii, which are described 
by the correlations. To see this, we use the small-/c behav- 
ior of the structure, S{k,t) = {tpkip-k) ~ R'^+^k'^ ©H- 
At time ti, the fluctuations SV^ in the volume covered 
by -0 = 1 within an region of size R2 is (up to numeri- 
cal factors) SV;^ - / d'^k J^' d\ J^^ d'^r' e' S{k, h) - 

R¥ lo^"" d'kSik,h), thus 6V^ ~ ±R['+^^/'Ri'-^^/'. 
In the case that the drops at t2 coincide with positive 
fluctuations in the drop density at ti, then this volume 
SV[p will be contributed to A{ti,t2) for every drop at 
t2, of number density l/R^- This gives the autocorre- 
lation decay A(ti,t2) (i?i/i?2)^'*^'*^^^, which saturates 
the YRD bound ||l|l. In the case of weaker correlations 
between drops at t2 and fluctuations in drop density at 
ti, the autocorrelations decay faster. With both cases, 
we recover the YRD bound A > d/2 + 2. Comparison 
with the e ^ 0+ result A = rf, we see that, at least for 
d < 4, the asymptotic decay of autocorrelations becomes 
faster when correlations suppress fluctuations in the drop 
density at large length-scales. 

In summary, we have demonstrated the existence of 
power laws in the autocorrelations, persistence and poi- 
soning of LSW and GC systems for T < Tc coarsening, 
and have calculated the exact asymptotic amplitudes and 
exponents. Our results are universal for isotropic systems 
such as polymer blends (before hydrodynamic regimes). 
While the exponent 9 is small in the dilute regime, it 
should be measurable since (|3|) holds for all t2 > ti. 

In the future, we hope that the poisoning exponents 
6*2 are measured for foams and large-g Potts models as a 
sensitive probe of their equivalence, and compared with 
(^. It would also be worthwhile to extend the indepen- 
dent interval approximation Q to the dilute globally- 
conserved case, to further explore the common 7^ ~ ^/d 
asymptote. Finally, our late-time crossover and logarith- 
mic corrections in the autocorrelations may survive in 
larger filling fractions. Studying A and its crossover as a 
function of e may resolve the current ambiguities [ p^pT| 
about autocorrelation decay in conserved systems. 
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